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Abstract 

In a previous paper, we reported the results of a partial-wave analysis of the pion- 
nucleon (7 tN) differential cross sections (DCSs) of the CHAOS Collaboration and 
came to the conclusion that the angular distribution of their -K + p data sets is incom¬ 
patible with the rest of the modern (meson-factory) database. The present work, 
re-addressing this issue, has been instigated by a number of recent improvements in 
our analysis, namely regarding the inclusion of the theoretical uncertainties when 
investigating the reproduction of experimental data sets on the basis of a given 
‘theoretical’ solution, modifications in the parameterisation of the form factors of 
the proton and of the pion entering the electromagnetic part of the irN amplitude, 
and the inclusion of the effects of the variation of the a-meson mass when fitting 
the ETH model of the nN interaction to the experimental data. The new analysis 
of the CHAOS DCSs confirms our earlier conclusions and casts doubt on the value 
for the 7 tN £ term, which Stahov, Clement, and Wagner have extracted from these 
data. 

PACS: 13.75.Gx; 25.80.Dj; 11.30.-j 


Key words: nN elastic scattering; nN £ term 


1 Introduction 


Stahov, Clement, and Wagner [1] have evaluated the pion-nucleon (ttN) £ 
term from the 7r ± p differential cross sections (DCSs) of the CHAOS Collab- 




oration [2,3]; the extracted value was 59 ± 12 MeV. In Ref. [4], which was 
available online (Nucl. Phys. A web site) a few months prior to the publica¬ 
tion of their paper [1], we had reported the details of a partial-wave analysis 
(PWA) of the same data (DENZ04, in our notation). To avoid influence from 
extraneous sources, an exclusive analysis of the DENZ04 DCSs had been per¬ 
formed in Ref. [4], after applying to these data the same analysis criteria 
which had been used earlier [5] in a PWA of the rest of the low-energy (pion 
laboratory kinetic energy T < 100 MeV) n ± p clastic-scattering measurements. 

In the Erst part of our paper [4] (optimisation), we had analysed the DENZ04 
DCSs following two theoretical approaches, one featuring standard low-energy 
parameterisations of the s- and p -wave AMnatrix elements, the other using the 
s- and p -wave A'-matrix elements of the ETH model 0 [6]. Both ways failed 
to produce acceptable results. We subsequently investigated the reproduction 
of the DENZ04 DCSs on the basis of the results of Ref. [5]. The compari¬ 
son revealed large discrepancies in the DENZ04 n + p data sets at forward and 
medium scattering angles, at all five energies covered by the CHAOS exper¬ 
iment. We thus concluded that the angular distribution of the DENZ04 ir + p 
DCSs was incompatible (in shape) with the rest of the modern (meson-factory) 
low-energy database. 

Owing to the fact that the present paper reports the results of a new analysis of 
the DENZ04 DCSs, the overlap with the material of Ref. [4] is inevitable. There 
are three main reasons for embarking on a new analysis of these measurements. 

• Included now in our results, for the first time, are also the uncertainties of 
the theoretical values (Syfj) when investigating the reproduction of experi¬ 
mental data sets on the basis of a given ‘theoretical’ solution (from now on, 
‘baseline solution’ or BLS); these uncertainties are sizeable in some kine- 
matical regions covered by the CHAOS experiment, e.g., in backward n~p 
elastic scattering. It is worth noting that, compared to our past PWAs, the 
uncertainties SyF are now somewhat larger, as they also contain the effects 
of the variation of the u-meson mass m a (see below). 

• Recent developments regarding the proton electromagnetic form factors sug¬ 
gest the replacement of the forms we have used earlier. The parameterisation 
of the Dirac Ff(t) and the Pauli Af (t) form factors of the proton with (tra¬ 
ditional) dipole forms has been found to provide a poor description of the 
‘world’ electron-proton (ep) unpolarised and polarised data [7]. Although 
the sensitivity of our results to the details of the parameterisation of these 
quantities is low (due to the smallness of the Q 2 transfer at low energies), we 


1 The ETH model of the irN interaction contains f-channel a- and p-exchange 
graphs, as well as the s- and u-channel contributions with all the well-established s 
and p baryon states with masses below 2 GeV; the model amplitudes obey crossing 
symmetry and isospin invariance. 
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will nevertheless adopt a more recent parameterisation scheme. In Ref. [8], 
the authors made use of the so-called Pade parameterisation [9] for the 
Sachs electromagnetic form factors G p E [t ) and G p M (t ) (in Ref. [8], the su¬ 
perscript p is omitted), and obtained the optimal values for the relevant 
parameters from a fit to ep measurements (see their Table II). The value of 
the proton rms electric charge radius >, evaluated from the content 

of that table (in fact, from the values of the parameters g 2 and g 6 ), is almost 
identical to the result of Ref. [10], thus disagreeing with the ^< r/ > value 
extracted from muonic hydrogen [11]. From now on, we will use the results 
of Ref. [8] in our analysis software. The pion form factor F n (t) is usually 
parameterised via a monopole form, e.g., see Ref. [12]. Despite the fact that, 
in the low-energy region, results of similar quality are obtained with either 
a monopole or a dipole form, we will adopt the monopole parameterisation 
henceforth. 

• The fits of the ETH model to the experimental data now involve the vari¬ 
ation of m a within the interval which is currently recommended by the 
Particle-Data Group (PDG) [13], namely between 400 and 550 MeV. The 
earlier fits were made assuming a fixed m a value (860 MeV [14]). 


There are three additional reasons, albeit less important, for revisiting this 
subject, a) The s-charmel contribution of the Roper resonance to the K -matrix 
element K 1 _ , as given in Section 3.5.1 of Ref. [15], is now explicitly included 
in the Erst step of the optimisation (low-energy parameterisations of the s- 
and p-wave K -matrix elements); this change induces very small effects in the 
fits to the 7 i~p elastic-scattering data, b) An improved approach for fixing the 
(small) d and / waves has been implemented; to suppress artefacts which are 
due to the truncation of small values, simple polynomials are now fitted to the 
d- and /-wave phase shifts of the current solution of the SAID analysis [16]. 
c) Used in Refs. [4,5] were the results of an earlier compilation of the physical 
constants by the PDG; the results of the most recent compilation [13] are used 
herein. 


The present paper is organised as follows. In Section 2, we give the details 
regarding the assessment of the goodness of the reproduction of an experi¬ 
mental data set on the basis of a BLS. The tests are now optimally structured 
and may be generally used in such investigations; for instance, we can use 
the prescribed tests to assess the reproduction of the n~p charge-exchange 
( 7 r~p —y 7 T°n) data on the basis of elastic-scattering results (pursuing the in¬ 
vestigation of the violation of the isospin invariance in the hadronic part of the 
nN interaction). Section 3 contains the results we obtained from the new anal¬ 
ysis of the DENZ04 DCSs. We briefly discuss the implications of our findings 
in Section 4. 
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2 Assessment of the goodness of the reproduction of an experi¬ 
mental data set on the basis of a BLS 


Prior to advancing to the technical details, we will give a few definitions which 
are of relevance in the discussion; if introduced later, they might disrupt the 
smooth description of the three tests delineated in the present section. For our 
purposes, a BLS is defined as follows. 

• A BLS is a set of values and associated uncertainties (yj^, Syfj , i G [1, Nj]), 
corresponding to the values of the kinematical variables, i.e., of the centre- 
of-mass scattering angle 6 and of T, at which the experimental data 
8i/ij P , i G [1, Nj]) have been acquired. The indices i and j identify the 
particular measurement, namely as the i th data point of the j th data set. 
The number of data points of the j th data set is denoted as N y 

• A BLS comprises predictions obtained via a Monte-Carlo simulation taking 
into account the results of the optimisation (i.e., the fitted values and the 
uncertainties of the model parameters, as well as the covariance matrix of 
each fit) of a PWA of nN data. 

Being a sum of independent normalised residuals, each following the normal 
distribution, our test-statistic is expected to follow the \ 2 distribution. As we 
will concentrate on this distribution hereafter, aiming at the identification of 
data sets which are poorly reproduced, we will tailor all expressions of the 
present section to one-sided tests (right-tail events). 

Let us assume that the background process, underlying the phenomenon under 
investigation, is a stochastic one (or, equivalently, that the null hypothesis 
is valid), described by the probability density function f(x) > 0, where x G 
[0, oo) is a numerical result obtained via a measurement made on the observed 
system. Kolmogorov’s second axiom dictates that 

p OO 

/ f(x)dx = 1 . (1) 

Jo 

The so-called p-value [f] is defined as the upper tail of the corresponding cu¬ 
mulative distribution function: 

poo 

P(Ao) = / f(x)dx . (2) 

Jxo 

The p-value is the probability that a measurement, obtained from the observed 
system, yield a result x which is more statistically significant than xq (in our 
case, that x > x 0 ). Assuming the validity of the null hypothesis, the p-value 


2 It is casual to refer to p-values in the statistical hypothesis testing in most domains 
of basic or applied research in economics, psychology, biology, medical physics, etc. 
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may thus be used as a measure of the result xq being due to chance: ‘small’ 
p-values attest to the statistical significance of that measurement. 

Of course, before assessing the statistical significance of a measurement, one 
must define what is meant by ‘small’ p-values. Unfortunately, the threshold 
signifying the outset of statistical significance is subjective; in reality, the 
setting of the significance level GO p mm rests on a delicate trade-off between 
two risks: a) of accepting the alternative hypothesis (of an effect not being due 
to statistical contrivance) when it is false and b) of rejecting the alternative 
hypothesis when it is true. Of relevance in the choice of the p min value is which 
of these two risks is being assigned greater importance. For instance, if the 
implications of risk (b) are considered to be more severe, compared to those 
of risk (a), an increase of the p m j n value is tenable. 

Most statisticians accept p min = 10~ 2 as the outset of statistical significance 
and p m .in = 5 • 10~ 2 as the threshold indicating probable significance. An in¬ 
teresting recent paper [17] interprets the lack of reproducibility of scientific 
results in various disciplines as evidence that the currently-accepted p min val¬ 
ues are rather ‘optimistic’; the author thus recommends the reduction of these 
thresholds by one order of magnitude P~l . 

The probability density function of the y 2 distribution with v > 0 degrees of 
freedom (DOF) is 


f(x,u) 


2 "l 4 {v/ 2 ) xVl2 1 exp(—x/2), for x > 0 
0 , otherwise 


( 3 ) 


where T(y) is the standard gamma function: 

poo 

r (y)= f 2 ' -1 exp(-f)dt . (4) 

Jo 

For a quantity x following the y 2 distribution, the expectation value E[x\ is 
simply equal to v and the variance E[x 2 ] — ( E[x ]) 2 is equal to 2v. The relation 
E[x\ = v has led many physicists to the use of the reduced y 2 value (i.e., of the 
ratio y 2 /z/) when assessing the goodness of the data description in modelling; 
as long as y 2 /v ~ 1, the results are claimed to be satisfactory. Of course, the 
interesting question in the statistical hypothesis testing relates to the value of 
y 2 /v at which the results start appearing unsatisfactory; evidently, a threshold 

3 The significance level is usually denoted as a in Statistics. 

4 Although Ref. [17] states that “nonreproducibility in scientific studies can be 
attributed to a number of factors, including poor research designs, flawed statistical 
analyses, and scientific misconduct”, we believe that, at least as far as the research 
in 7 tN physics is concerned, the main reason might simply be ‘excessive optimism’ 
when assessing the systematic effects in the experiments; in all probability, these 
uncertainties are frequently, if not systematically, underestimated. 
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value for y 2 t v ma y be extracted from p min , yet it turns out to be //-dependent, 
hence cumbersome to use. Such a departure from simplicity is meaningless. It 
makes more sense to perform the direct test and assess the statistical signifi¬ 
cance by comparing the p-value, associated with the observed y 2 f° r v DOF, 
with p mint this is achieved by simply inserting f(x, v) of Eq. (3) into Eq. (2), 
along with xo = y 2 , and evaluating the integral; several software implemen¬ 
tations of dedicated algorithms are available, e.g., see Refs. [18] (Chapter on 
‘Gamma Function and Related Functions’) and [19], the routine PROB of the 
CERN software library, the functions CHIDIST/CHISQ.DIST.RT of Microsoft 
Excel, etc. 

One additional ‘side’ remark is due. Various definitions of the ‘data set’ have 
been in use, involving different choices of the experimental conditions which 
must remain stable/constant during the data-acquisition session. The proper¬ 
ties of the incident beam, as well as the (physical, geometrical) properties of 
the target, have been used in the past in order to distinguish the data sets 
of experiments performed at one place (i.e., at a meson factory) over a short 
period of time (typically, a few weeks). However, data sets have appeared 
in experimental reports relevant to the nN system, which not only involved 
different beam energies, but also contained measurements of different reac¬ 
tions (e.g., mixing n + p and tt~ p elastic-scattering measurements). As a result, 
the only prerequisite for accepting measurements as comprising one data set 
is that they share the same absolute normalisation (and, consequently, nor¬ 
malisation uncertainty Szj). Of course, this is only a prerequisite, hence a 
necessary, not a sufficient, condition. The decision regarding the acceptance 
of a set of measurements as comprising one data set cannot be made without 
an investigation of the stability of the experimental conditions at which the 
raw measurements had been acquired (this may be difficult to assess), as well 
as of their (on-line and off-line) processing on the way to the extraction of the 
final experimental results. 

We now enter the details of the reproduction. Let us assume that the absolute 
normalisation of the j th data set is known up to a relative uncertainty Szj. 
Let us also assume that none of the important quantities, appearing in the 
denominators of the expressions of the present section, vanishes. (For the sake 
of compatibility with our past works [4,5,20], we will retain the index j in the 
expressions, despite the fact that its use in this section is redundant.) 

One way of assessing the goodness of the reproduction of the j th data set 
involves the determination of the amount of scaling (application of a multi¬ 
plicative factor to the BLS, enabling the increase or decrease of its values, 
resulting in its ‘upward’ or ‘downward’ shift as ‘one piece’) which must be ap¬ 
plied to the BLS (yjj, Syf-, i G [1, At,-]) in order that it ‘best’ accounts for the 
entire data set (y®J p , Sy^j P , i G [1, N 3 ]). Regarding the reproduction of data 
sets by a BLS, we have gained inspiration from the Arndt-Roper formula [21], 
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which we have been using in our fits to the experimental data since Ref. [20]. 
We will next propose tests of the overall reproduction, of the shape, and of 
the absolute normalisation of the data sets. 


We first evaluate the ratios, r l3 = y^ p /yjj] if the quantities y f - 3 P and yjj are 
independent (which is certainly true in our case, as the DENZ04 measurements 
have not been used in the determination of the BLS), the uncertainties hry,- 
are obtained via the application of Gauss’s error-propagation formula: 


Snj = fij 


\ 


exp 


2 

+ 


/MV 

V vij ) 


( 5 ) 


The goodness of the reproduction is assessed on the basis of the function Xj( z j) 
defined as: 


n. 


Xj(zj) 


= E 

i=l 


>'ij -/ 
v 5r a , 


+ 


z i - 1 

v 5z j , 


( 6 ) 


It is convenient to introduce the weights Wij via the relation Wij = (hr^) 2 


The second term on the right-hand side of Eq. (6) takes account of the ap¬ 
plication of scaling to the BLS. This contribution depends on how well the 
absolute normalisation of the j th data set is known. For a poorly-known abso¬ 
lute normalisation, Szj is large and the resulting contribution from the scaling 
is small; the opposite is true in case of a well-known absolute normalisation. 
Evidently, the ‘best’ reproduction of the j th data set is achieved when, by 
varying the scale factor Zj, the function xj( z j) is minimised: 


MM = 0 

dzj 


( 7 ) 


The solution of this equation is 


>7 V 'i iriji-ij - (dr:/) 2 


( 8 ) 


Inserting this expression for Zj into Eq. (6), one obtains 


1 Nj Nj Nj 2 

(Xj)min = WE 7Tx wf ( W V V W V r ij ~ ( X] 

Ei=i + {ozj) - v i=1 i=1 i=1 

Nj 

■b {$ z j) Wj 3 (rn — i) 

i= 1 



( 9 ) 


Expression (9) yields the minimal y 2 value in the description of the j th data 
set, containing Nj data points. In fact, one additional measurement had been 
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made on this data set, namely the one fixing its absolute normalisation, which 
is known with relative uncertainty Sz r As a result, the number of DOF for this 
data set is equal to Nj + 1 — 1 = At,-; the subtraction of one unit is due to the 
use of Eq. (8) as a constraint, fixing the value of the scale factor Zj. Therefore, 
the quantity (y 2 ) m m of Eq. (9) is expected to follow the y 2 distribution with 
v = Nj DOF. To obtain the p-value of the overall reproduction of the j th data 
set, one uses Eq. (2) with f(x ) = f(x, u) of Eq. (3), along with Xo = (Xj)min 
and v = Nj. 

Two additional tests on each data set are possible. These tests are particularly 
useful in case that the overall reproduction of a data set is poor; they point 
to the culprit for the poor overall reproduction, namely to the shape or to the 
absolute normalisation of the data set. 

• To examine the shape of the j th data set (with respect to that of the BLS), 
it is important to allow the BLS to reproduce the data set ‘optimally’, i.e., 
regardless of the scaling. This is equivalent to setting Szj —* oo or (Szj)~ 2 = 
0 in Eqs. (8) and (9). The corresponding quantities will be denoted as Zj 
and (Xj)stat, respectively; the quantity (Xj)stat represents the fluctuation in 
the j th data set which (assuming the correctness of the shape of the data 
set) is of pure statistical nature. 


z i = 


^Nj 

11; I ll'ijrij 

^Nj 

Eiii Wij 


(Xj)stat ^Ni 


J2 "A X w a r ij 


Hi —1 IVij 8 j = ] i=l 



As expected, both expressions are identical to those derived for the weighted 
average of a set of independent measurements, as well as for the correspond¬ 
ing y 2 value for constancy. Owing to the fact that the normalisation uncer¬ 
tainty is not used in Eq. (11), the quantity (xf)stat is expected to follow the 
y 2 distribution with u = Nj — 1 DOF. The p-value, obtained from Eq. (2) 
with x 0 = (Xj)stat an d v = Nj ** 1, may be used in order to assess the 
constancy of the input values ry,- or, equivalently in our case, to examine 
the shape of the j th data set with respect to the BLS 0- 
• To assess the compatibility of the absolute normalisations of the j th data 
set and of the BLS, one first estimates the scaling contribution to (y 2 ) r „m 


5 In fact, the test simply assesses the goodness of the representation of the input 
data by one overall average value. A failure indicates either a bad shape (e.g., a slope 
being present in the input data) or ‘scattered’ input values with small uncertainties. 
Visual inspection of the input data reveals that the latter option is not the case in 
the problem we have set out to investigate. 





via the relation: 




min 



(Szj) 2 I 

(e£i Wijir-i 

> -1)) 

2 

(e£i<a 

ij + (Vi)- 2 ) 

i Wij 


( 12 ) 


where use of Eqs. (9) and (11) has been made. The quantity (y 2 ) sc is ex¬ 
pected to follow the x 2 distribution with 1 DOF (which, of course, is the 
normal distribution). 


We now summarise the three tests of the goodness of the reproduction of the 
j th data set by a BLS. 

• The overall reproduction can be tested using (Xj)min °f Eq. (9) as x 0 in 
Eq. (2) and v = Nj DOF. 

• The shape (statistical fluctuation) can be tested using (x 2 )stat of Eq. (11) 
as Xq in Eq. (2) and v = Nj — 1 DOF. 

• The absolute normalisation can be tested using (y 2 ) sc of Eq. (12) as x 0 in 
Eq. (2) and v — 1 DOF. 

The tests outlined in the present section are objective. The only subjective as¬ 
pect in the analysis pertains to the choice of the p min level signifying statistical 
significance. 


3 Results of the new analysis of the CHAOS DCSs 


We now present the results of the new analysis of the DENZ04 DCSs, using the 
changes in our approach as detailed in Sections 1 and 2. We first investigated 
the description of the experimental data on the basis of the standard K -matrix 
low-energy parameterisations which we employ as first step in our PWAs, 
identifying and removing any outliers from the database; in this step of the 
optimisation, the theoretical constraint of crossing symmetry is not imposed 
onto the fitted scattering amplitudes. We used the split data sets, as defined 
in Refs. [4,16]. For the truncated 7 r + p database, the minimal y 2 value was 
397.4 (to be compared to 401.2 in Ref. [4]) for 260 DOF. Exempting two data 
points which changed status @, the list of outliers is the same as in Ref. [4], The 
isospin-1 amplitudes were fixed from the final fit to the truncated n + p database 
and were imported into the analysis of the it~ p elastic-scattering database; 
only three data points were removed, the same measurements which had to 

6 The 19.90 MeV, 42.75° DCS, which was an outlier in Ref. [4], became an accepted 
data point. On the contrary, the 32.00 MeV, 40.83° DCS, which was an accepted 
measurement earlier, turned into an outlier. Both points are close to the acceptance 
threshold and even small changes in the approach may affect their status. 
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be excluded in Ref. [4]. For the truncated tt~ p elastic-scattering database, the 
minimal y 2 value was 352.9 (to be compared to 350.2 in Ref. [4]) for 261 DOF. 
Similarly to Ref. [4], the result for a cc , obtained in the final fit to the tt~ p 
elastic-scattering database, was unacceptable, about 0.042 m” 1 (m c denotes 
the mass of the charged pion), differing from the experimental result obtained 
directly at the ttN threshold [22,23] by a factor of 2; this discrepancy is dne 
to the inadequacy of the isospin-f amplitudes to account simultaneously for 
both elastic-scattering reactions. The common fit to the truncated combined 
7 r ± p elastic-scattering databases yielded no further outliers and a minimal y 2 
value of 750.4 (to be compared to 751.5 in Ref. [4]) for 521 DOF; the final a cc 
value was almost identical to the result we had obtained in the previous step, 
from the fit to the truncated tt~ p elastic-scattering database. 

The fits of the ETH model to the truncated combined n ± p elastic-scattering 
databases of the CHAOS Collaboration were next attempted, as described 
in Ref. [15]. The a-meson mass m a was varied within the interval which is 
currently recommended by the PDG [13]. Unfortunately, none of these seven 
(see Ref. [15], p. 178) fits terminated successfully. In all cases, negative diagonal 
elements were detected in the covariance matrix, urging the MINUIT software 
library [24] (FORTRAN version) to enforce positivity by adding arbitrary 
constants to the diagonal. As a result, the fitted values and the uncertainties 
of the model parameters are meaningless in all seven attempts to account for 
the DENZ04 DCSs on the basis of the ETH model. 

We also followed the recommendation of the CHAOS Collaboration and anal¬ 
ysed their unsplit (original) data sets. As explained in Ref. [4], one expects 
that the problems, which we encountered in the analysis of their split data, can 
only be aggravated when using the unsplit data sets. To start with, 72 data 
points were identified as outliers and had to be removed from the database: 
65 of these measurements belong to the n + p data sets, 7 to the n~p clastic- 
scattering ones. The result for a cc , obtained in the final fit to the resulting 
truncated combined tt ± p elastic-scattering databases using our low-energy pa- 
rameterisations of the s- and p- wave 77-matrix elements, was equal to about 
0.056 m~ l . We subsequently attempted to fit the ETH model to the truncated 
combined n ± p elastic-scattering databases. As in the case of the split data, all 
seven fits failed. Adding to the severity of the problems, as reported earlier 
when using the split data as input, the coupling constant g W NN came out in 
the vicinity of 0 in all seven attempts. It is not possible to obtain anything 
reasonable from the unsplit data sets of the CHAOS Collaboration. 

Given the seriousness of the problems we have encountered in analysing the 
DENZ04 DCSs, we can only investigate their reproduction on the basis of ex¬ 
isting BLSs. The BLS in the present paper is obtained from the results of new 
fits to the truncated combined n ± p elastic-scattering databases of Ref. [15], 
after applying the changes as detailed in Section 1; the differences to the re- 
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suits of Ref. [15] were small, not exceeding 40% of the quoted uncertainties, 
save for the model parameter n p where the difference amounts to about one 
standard deviation. The predictions for the ‘theoretical’ values yjj and their 
uncertainties Syf-' were obtained on the basis of 14 million Monte-Carlo events 
for each data point in the CHAOS database. As always, the routine CORGEN 
of the CERN software library was used in the generation of the Monte-Carlo 
events; input to CORGEN is the ‘square root’ of the covariance matrix, which 
was obtained from the optimisation results with the routine CORSET. The 
CPU consumption per energy and angle value was about 2 min on a fairly- 
fast personal computer. As they are obtained in an analysis of a larger number 
of measurements, the uncertainties 5yjj are generally expected to be smaller 
than the experimental uncertainties of single experiments. Nevertheless, they 
are sizeable in some kinematical regions covered by the CHAOS experiment, 
e.g., in backward n~p clastic scattering (see caption of Fig. 2 of Ref. [4]). 

We now present the results of the analysis of the ratios r t] = y^ p /yjj for each 
data point in the DENZ04 database. Ideally, these ratios are constant (i.e., 
independent of 9 and of T) and equal to 1 (reflecting perfect agreement in the 
absolute normalisation of the data set being tested with respect to the BLS). 
Constant values of r^, different from 1, point to differences in the absolute 
normalisation of the two sets of values, whereas a statistically-signihcant de¬ 
parture of rij from constancy is evidence of a discrepancy in the shapes of 
the two angular distributions of the DCS. The ratios r l3 for the DENZ04 7r ± p 
elastic-scattering databases are shown in Figs. 1 and 2. 

The p-values of the reproduction of the DENZ04 DCSs [2,3] are given in Table 
1; the table corresponds to the unsplit data sets. The eleven outliers, which had 
been identified in the first step of the analysis (i.e., employing standard low- 
energy parameterisations of the s- and p- wave Jl-matrix elements), have not 
been removed. The removal of these measurements induces very small effects 
and cannot alter the conclusions in the case of the n + p data; as the free scaling 
(free floating) of the backward-angle 25.80 MeV data set had been suggested 
in the Erst step of the optimisation, the treatment of the unsplit n + p 25.80 
MeV data set is not straightforward. The removal of the two outliers from the 
7 r~p 25.80 MeV data set improves its reproduction by the BLS; the problem 
with this data set mainly rests with the peculiar shape of the forward-angle 
measurements (the same behaviour is observed in the corresponding n + p data 
set), where the dominant contribution to the nN scattering amplitude comes 
from the electromagnetic interaction. 

All six p-values for the overall reproduction of the n + p data sets in Table 1 are 
smaller than the significance level p mm for the rejection of the null hypothesis 
which we adopt in our PWAs, namely p m i n ~ 1.24 • 10~ 2 ; this value of p m * n , 
equivalent to a 2.5 o effect in the normal distribution, is close to 1% which 
most statisticians adopt in the statistical hypothesis testing. The 37.10 MeV 
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7 r + p data set is the best reproduced by the BLS (p ~ 3.65 • 10 -3 ), whereas the 
25.80 MeV ir + p data set is obviously the worst reproduced. Beyond doubt, the 
problems in the reproduction of these data sets relate to their shape. On the 
other hand, the only n~p r tJ ratios, which are poorly reproduced, are those of 
the 25.80 MeV data set, which appears to contain suspicious measurements at 
small 6] the reproduction of the remaining ir~p elastic-scattering data sets is 
satisfactory. Evidently, the angular distribution of the DENZ04 7i + p database 
disagrees (in shape) with the rest of the 7r + p database, whereas the DENZ04 
7 T~p elastic-scattering DCSs appear to be in reasonable agreement with the 
BLS. We had reached the same conclusion in Ref. [4], Regarding the split data, 
four (out of 17) 7 r + p data sets are poorly reproduced, as is the forward-angle 
25.80 MeV Ti~p elastic-scattering data set. Given that our results now contain 
all contributions to the uncertainties, there is hardly room for improvement 
in the reproduction of the DENZ04 ir + p database. 


4 Discussion and conclusions 


In the present paper, we analysed the differential cross sections (DCSs) of the 
CHAOS Collaboration [2,3] and investigated their reproduction on the basis of 
the results obtained from a partial-wave analysis of the rest of the low-energy 
(pion laboratory kinetic energy T < 100 MeV) n±p elastic-scattering data 
[15]. Since our previous analysis of the same data [4] appeared, there have 
been three main developments calling for a new investigation of this subject. 

• When investigating the reproduction of an experiment on the basis of a 
given ‘theoretical’ solution, our results also include now the uncertainties 
Sylj of the theoretical values (see Eq. (5)). 

• Recent developments regarding the proton electromagnetic form factors sug¬ 
gest the replacement of the forms we used in our earlier analyses; from now 
on, we will adopt the parameterisation (and the optimal parameter val¬ 
ues) of Ref. [8]. Additionally, the pion electromagnetic form factor will be 
parameterised via a monopole form. 

• The fits of the ETH model [6,15] to the experimental data now involve the 
variation of the cr-meson mass (m a ) within the interval which is recom¬ 
mended by the Particle-Data Group [13] (currently, between 400 and 550 
MeV). 

The new analysis of the CHAOS DCSs demonstrated that the conclusions of 
Ref. [4] hold and that the angular distribution of their ir + p cross sections is in 
conflict with the rest of the modern (meson-factory) low-energy n + p database. 

One might argue that a disagreement between any two sets of data attests only 
to the faultiness of at least one of them. In order to determine which of the two 


12 



sets of values is flawed, the insight gained from theory is often valuable. To this 
end, we analysed the two sets of measurements, irrespective of one another, 
by employing two theoretical approaches in the analysis, namely the standard 
low-energy parameterisations of the s- and p -wave K -matrix elements and the 
ETH model. We are critical of the DCSs of the CHAOS Collaboration because 
we have not been able to obtain anything reasonable from the analysis of 
these data following either theoretical approach. On the contrary, we did not 
encounter problems when analysing the rest of the modern database following 
both theoretical approaches. As the beam energy in the CHAOS experiments 
was sufficiently low, an investigation of the description of their DCSs within 
the framework of the Chiral Perturbation Theory (e.g., with the method of 
Ref. [25]) should be possible. This would be an interesting subject to pursue. 
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Table 1 


The details of the reproduction of the differential cross sections of the CHAOS 
Collaboration [2,3] by the baseline solution (BLS) obtained following the approach of 
Ref. [15]. The results now include the uncertainties dyf- of the BLS (see Eq. (5)). The 
columns represent: the pion laboratory kinetic energy T (in MeV), the number of 
data points Nj of the j th experimental data set, and the three p-values corresponding 
a) to the overall reproduction of the data set, b) to the reproduction of its shape, 
and c) to the reproduction of its absolute normalisation. The table corresponds to 
the unsplit (original) data sets of the CHAOS Collaboration. 


T 


Overall 

Shape 

Absolute normalisation 

7 r + p scattering 

19.90 

33 

2.37 • 10“ 9 

1.30 • 10" 9 

8.46 • 10' 1 

25.80 

43 

9.97 • 10~ 78 

3.28 • 10~ 78 

6.00 • 10” 1 

32.00 

46 

2.18 • 10“ 4 

1.53 • 10 -4 

9.67- 10” 1 

37.10 

49 

3.65 • 10“ 3 

4.99 • 10" 3 

1.00 • 10' 1 

43.30 

53 

3.19 • 10” 11 

1.87 • 10 -11 

8.65 • 10” 1 

43.30(rot.) 

51 

1.16 • 10~ 12 

7.22 • 10” 13 

5.81 • 10’ 1 

7 T p elastic scattering 

19.90 

31 

9.90 • nr 1 

9.89 • 10" 1 

4.38 • 10” 1 

25.80 

45 

4.82 • 10~ n 

3.76 • 10 -11 

3.30 • 10’ 1 

32.00 

45 

1.40 • 10“ 2 

5.26 • 10 -2 

4.31 • 10~ 3 

37.10 

50 

3.55 • 10 _1 

4.08 • 10" 1 

1.17 • 10” 1 

43.30 

51 

1.44 • nr 1 

2.47- 10" 1 

2.12 • 10“ 2 

43.30(rot.) 

49 

6.38 • 10’ 1 

8.71 • 10 _1 

5.27 • 10~ 3 
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Fig. 1. The -K + p differential cross sections of the CHAOS Collaboration [2,3] 
normalised to the corresponding predictions ( yjj ) obtained following the approach of 
Ref. [15]; the eight outliers in the ir + p database are also contained in the figure. The 
normalisation uncertainties of the experimental data sets (see Refs. [2,4] for details) 
are not shown. Unlike Ref. [4] (e.g., see caption of Fig. 1 therein), the uncertainties 
8yfj are now included in the results (see Eq. (5)). 
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Fig. 2. The n p elastic-scattering differential cross sections of the CHAOS Col¬ 
laboration [2,3] normalised to the corresponding predictions (yjj) obtained 

following the approach of Ref. [15]; the three outliers in the tt~ p elastic-scattering 
database are also contained in the figure. The normalisation uncertainties of the ex¬ 
perimental data sets (see Refs. [2,4] for details) are not shown. Unlike Ref. [4] (see 
caption of Fig. 2 therein), the uncertainties Syjj . which are sizeable at backward 
angles, are now included in the results (see Eq. (5)). 
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